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Recent angle-resolved photoemission and neutron scattering data have provided new ingredients 
for the interpretation of scanning tunneling spectra on F^S^CaC^Og+i- We analyze the low- 
temperature tunneling spectra, from oxygen overdoped to underdoped samples, including details 
about the bilayer splitting and the neutron resonance peak. Two van Hove singularities are identified: 
the first is integrated in the coherence peaks, the second is heavily broadened at higher binding 
energy. The shape of the tunneling spectra suggests a strong coupling of the quasiparticles with a 
collective mode, and a comparison with photoemission shows that the scattering rate in tunneling 
is an order of magnitude smaller than in ARPES. Finally, the theoretical spectra calculated with an 
isotropic tunneling matrix element are in better agreement with the experimental data than those 
obtained with anisotropic matrix elements. 

PACS numbers: 74.50.+r, 73.40.Gk, 74.72.-h, 74.72.Hs 



I. INTRODUCTION 

Scanning tunneling spectroscopy (STS) is a powerful 
tool to study the electronic properties of solids. Its re- 
markable energy and spatial resolution makes it particu- 
larly well suited for materials characterized by small en- 
ergy and short length scales, like the cuprate high-T c su- 
perconductors (HTS). Among the HTS, the bilayer com- 
pound Bi 2 Sr 2 CaCu 2 0s+,5 (BSCCO) has often been stud- 
ied by STS, because it cleaves easily and offers an atom- 
ically flat BiO surface. It is possible to tunnel through 
the insulating BiO and SrO surface layers into the Cu02 
plane, where all the exciting properties of the cuprates 
are believed to reside. Many important results on the 
nature of the superconducting, normal, ai ^WHWj4rfj a te 
of BSCCO have been obtained using STS.Mffiliai 

In the superconducting state, the main feature of the 
differential tunneling conductance (dl/dV) spectrum is 
the quasiparticle excitation gap, which has been observed 
in BSCCO stod-jSljudied as a function of doping and 
temperatureoB'BMl-j The presence of excitations within 
the superconducting gap, linearly increasing with energy 
around V = 0, indicates that the order parameter has 
nodes, and presumably d x 2_ y 2 symmetry. Other char- 
acteristics of the BSCCO spectra are the celebrated dip- 
hump structure at energies larger than the excitation gap, 
an asymmetry between electron and hole tunneling, and 
coherence peaks with considerably more spectral weight 
than predicted by BCS theory. None of these three char- 
acteristics has been met with an explanation that is com- 
monly agreed upon. 

Under certain assumptions, Bill and considering the 
nearly two-dimensional nature of BSCCO, one finds that 
the shape of the STS spectrum is determined by only 
three ingredients: the bare electron dispersion in the 
Cu0 2 plane, the self-energy S(fc,w) which embodies all 



electronic interactions, and a tunneling matrix element 
Tfc which couples the electronic states of momentum k 
at the sample surface with the metallic tip. In the par- 
ticular case where Tk is a constant, the STS spectrum di- 
rectly relates to the quasiparticle density of states (DOS) 
in the Cu0 2 plane. Behind this apparently simple state- 
ment lies a real difficulty to disentangle the contribu- 
tion of each of the above ingredients. Various features in 
the STS spectra were thus attribut ed e ither to van Hove 
singularities in the band-structure£-st3 or to self-energy 
effects.EErEZI Furthermore, the general shape of the low- 
temperature spectrum was. C | On .s ideij ed suggestive of an 



anisotropic matrix element .t3E3E3E3 

The situation has recently become even more compli- 
cated with the observation of a clear bila«eiusalitting by 
angle-resolved photoemission (ARPES) H'l3E3 The two 
Cu0 2 layers in the BSCCO unit cell give rise to two non- 
degenerate bands close to the Fermi energy. As a conse- 
quence, there are two bands — instead of one, as assumed 
previously — which can contribute to the dl/dV spectra. 
Using ARPES, the shape of these bands has been deter- 
mined in_the normal state of underdoped and overdoped 
BSCCO.Efl 

In fact, the correct interpretation of STS spectra relies 
upon a realistic modeling of the data, and different mod- 
els can lead to opposite conclusions. Moreover, the mocL. 
eling critically depends on details of the band structurco 
and of the spin excitation spectrum r^r^ which have be- 
come available only very recently. 

In this study we compare the predictions of several 
models to the STS spectra measured at low tempera- 
ture an BSCCO samples with different oxygen dopings 
levels.!!! Our calculations take into account the bilayer 
splitting and are based on the band structures deter- 
mined in Ref. ^4]. We assume a pure d x 2_ y 2 symmetry 
of the superconducting gap and consider three different 
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models for the self-energy: a conventional BCS model; 
a pheaomenological model proposed to fit the ARPES 
data,E-3 thus allowing a direct comparison between pho- 
toemission and tunneling; and a model which describes 
the coupling of quasiparticles to a collective mode. Ell We 
also compare the effects of isotropic and anisotropic tun- 
neling matrix elements. 



II. MODELING THE TUNNELING SPECTRA 

In the tunneling-Hamiltonian formalism,^ the sample 
and tip are coupled by a matrix element Mk q which rep- 
resents the overlap of the electronic states on both sides 
of the tunnel junction. The resulting differential conduc- 
tance at bias voltage V is given by 



dl_ 

dV 



k.n 



T k \ 2 A n (k,u)f{Lu 



eV) (1) 



where / is the Fermi function and A n is the spectral func- 
tion in the sample. The sign convention is such that, at 
zero temperature, negative energies to and eV correspond 
to occupied states, and positive energies to unoccupied 
states. The new matrix element appearing in Eq. (|l|) is 
\Tk\ 2 = Y^q \Mkq\ 2 A t i p (q,uj). Using a tip with a feature- 
less DOS, we can assume that it is energy independent. 
According to Ref. [IJ the dependence of Tk on k z is can- 
celled by the band dispersion along k z [z is the tunneling 
direction). In Eq. (|l|), n refers to the two bands result- 
ing from the bilayer splitting as discussed below. We 
assume that the k xy dependence of Tk is the same for 
both bands; this then leads to an n-independent matrix 
element. A n is related to the electron dispersion and 
self-energy through 



A(k,oj) = --Im 

7T 



1 



£r-e fe -£(fc,cj) 



(2) 



The band index is omitted for simplicity. The lifetime 
broadening T is introduced here for computational con- 
venience and is set to T = 1 meV in all of our calculations. 

To account for the bilayer splitting, we use the band 
structures determined in Ref. for the anti-bonding (A, 
plus sign) and bonding (B, minus sign) bands: 



-(A,B) 



—2 t(cos k x + cos k y ) + 4 1' cos k x cos k y 
—2 t"(cos 2k x + cos 2k y ) 
±jt^(cos k x — cosky) 2 + Ae. 



(3) 



The interlayer coupling is described by t± ; the maximum 
energy splitting between the A and B bands is 2t± and 
coincides with the van Hove singularities at the (tt, 0) 
point in the Brillouin zone. The tight-binding parame- 
ters inferred from ARPES for overdoped (OD69.0K) and 
underdoped (UD77.0K) BSCCO are given in Table Q. 
These parameters deviate slightly from those reported in 
Ref. pi t± was determined in Ref. M for the OD69.0K 



sample. The determination of t± for the UD77.0K sample 
is complicated by the strong influence of the pseudo- and 
superconducting gaps on the photoemission spectra. For 
that reason, UD77.0K was assumed to have the same t± 
as OD69.0K in Ref. ^4|. In this work, we have estimated 
t± for UD77.0K from a careful comparison between lead- 
ing edge gaps of the A and B bands. E3 The remaining 
parameters were modified accordingly. This leads to a 
better agreement with the tunneling data. From the pa- 
rameters of OD69.0K and UD77.0K, we make a linear 
interpolation with doping to obtain the parameters-ap- 
propriate for the samples studied by STS (Table |) xB 

The first model we consider is the conventional BCS 
model. In this case the self-energy reads: 



S (1) (fc,w) = 



I A* 



iT + e k ' 



(4) 



where the gap has d x i_ y i symmetry: A*. = ^ A (cos k x — 
cosky). The only free parameter, Ao, is adjusted to the 
experimental data. 

In the second modaL we follow an attempt to fit 
ARPES measurementscil by complementing the BCS 
model with a phenomenological fc-independent self- 
energy T,q(lu). At low energy, So(w) has a marginal Fermi 
liquid form: 



EoM = -Aco - W{a u) 2 + (Pok B T) 2 . 

Good agreement with ARPES was found— with, 
rameters A = 1, a = 2, and A) = 



if 



(5) 

pa- 
The 



complete self-energy including the superconducting gap 
is 



£( 2 >(fc,^) = I]o(a 



|A fc | 2 



OJ + iT + Bk - So(w) 



(6) 



Eq. ^ can be recast in a form similar to Eq. (^), with 
the bare dispersion Ek and gap A& renormalized to cek 
and cAfc, respectively, and the lifetime broadening Y re- 
placed by r — clmEo(w) with c = (1 + A) -1 . We will use 
Eq. @ to check if the self-energy inferred from ARPES 
is compatible with the STS data. 

Finally, the third model we consider accounts for the 
possible interaction of the quasiparticles with a collective 



TABLE I: Tight-binding parameters for the conduction bands 
of BSCCO. The parameters for OD69.0K and UD77.0K follow 
from fits to ARPES data, and in the other cases from a linear 
interpolation with doping. All numbers are in eV. 



Sample 


t 


t' 


t" 


tx 


Ae 


OD69.0K 


0.40 


0.090 


0.045 


0.082 


0.431 


UD77.0K 


0.39 


0.076 


0.034 


0.097 


0.304 


OD56.0K 


0.40 


0.092 


0.047 


0.081 


0.449 


OD74.3K 


0.40 


0.089 


0.044 


0.084 


0.422 


OP92.2K 


0.39 


0.082 


0.039 


0.090 


0.361 


UD83.0K 


0.39 


0.077 


0.035 


0.096 


0.317 
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mode. Such a mode was indeed observed in BSCCO bjj 
neutron scattering, with momentum centered at (71", ir).c3 
It appears at an energy Q « 5.4 ksT c f3 and is character- 
ized by a correlation length £ < 2 in lattice units (we will 
take £ = 2 in our simulations). The interaction of this 
mode with the electrons was described theoretically in 
Ref. and involves a coupling constant g for which we 
choose the value g = 0.65 eV. The self-energy entering 
Eq. (0) is 



OD56.0K 
A =20meV 
o exp. data 



E (3) (fc,w) = E u (fc,w) + 



|A fc + £ 12 (fc,o/)| 2 

- «r + £k - £22 (fc, w ) 



(7) 



The components E« are a convolution of the bare BCS 
propagator (in Nambu representation) with the spin sus- 
ceptibility The latter is represented by a simple analyt- 
ical function which approximates the neutron measure- 
ments. We refer the reader to Ref. ^Tjfor further details. 

Very little is known about the actual fe-dependence 
of the tunneling matrix element, although band calcu- 
lations suggest that it is anisotropic with a shape ccp, 
responding to the dispersion of the bilayer splitting,L3 
i.e. (cos k x — cos k y ) 2 . Such a matrix element would pro- 
hibit tunneling into the nodal direction (jr, tt) and would 
highlight the region (71", 0) of the van Hove singularities. 
We will consider the two limiting cases of a completely 
isotropic T s = Tq and anisotropic Td = To (cos k x — cos k y ) 
matrix element, as well as admixtures of the form |Ifc| 2 = 
a\T s \ 2 + (1 - a)\T d \\ 

For each model we calculate the differential tunneling 
conductance using Eq. (0) and a 1024 x 1024 mesh of k 
points. The temperature is set to 4.2 K. The experimen- 
tal energy resolution is simulated by filtering the data 
with a Gaussian of width 1 meV. Finally, unless stated 
otherwise, the calculated spectra are normalized to the 
total spectral weight of the experimental data over the 
energy range from —300 to +300 meV. 



III. RESULTS AND DISCUSSION 

A. BCS model 

To compare the calculated spectra with experimental 
results over a wide doping range, we use low-temperature 
(4.2 K) data of Renner et alM The results of the BCS 
model, with an isotropic matrix element, are shown in 
Fig. |l|. The agreement at subgap energies is good, except 
for OD74.3K where the spectral weight is somewhat un- 
derestimated. The (calculated and experimental) spec- 
tra show the V-shape at zero bias typical for a e?-wave 
order parameter. At higher energies, however, the model 
misses the dip-hump feature, which is most pronounced 
in the optimally doped and underdoped spectra at —60 to 
— 100 meV, and the model generally presents too sharp 
structures. In particular, the van Hove singularities of 
the A and B bands show up unrenormalized in the spec- 
tra while no such singularities exist in the experimental 




-300 -200 -100 100 200 300 
V (mV) 

FIG. 1: Tunneling spectra in the conventional BCS model, 
with a d-wave superconducting gap and an isotropic matrix 
element. The contributions of the A and B bands are shown 
separately, and their sum should be compared to the experi- 
mental data (circles). The arrows roughly indicate the maxi- 
mum of the broad feature in the background (see text). 



curves. None of these discrepancies can be reduced by the 
inclusion of an anisotropic matrix element, which would 
suppress spectral weight below the gap energy and raise 
the van Hove singularities even more. Nevertheless, with 
this first approach we can draw some conclusions which 
remain valid for the more sophisticated models discussed 
below. 

From the difference between the A and B bands, it 
becomes clear that most of the weight of the coherence 
peaks in the tunneling spectra is related to the van Hove 
singularity of the A band, which lies close to the Fermi 
level. The presence of the van Hove singularity, in the 
coherence peaks explains their unusual height Jia Note, 
however, that without the B band the peaks would be- 
come too high (with respect to the background) to fit 
the experimental data. The van Hove singularity also 
explains why the peak at negative bias is generally a bit 
higher than the one at positive bias. 

In general, we observe that the van Hove singularity of 
the A band remains close to the Fermi level (integrated 
in the coherence peaks) for all doping levels. The van 
Hove singularity of the B band, however, moves to higher 
binding energy with underdoping. This shift is related to 
the combined effects of increasing gap and increasing ij_ 
with underdoping. It is in qualitative agreement with 
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the behavior-observed in the background of the experi- 
mental datafra on going from overdoped to underdoped, 
the background of the spectra becomes more asymmetric. 
One can distinguish a very broad feature moving away 
from the Fermi level with underdoping. In fact, we have 
shown that this background can be fitted very well by 
a broadened van Hove singularity at energies consistent 
with the doping level, in a rigid single-band picture.E£l 
These energies are indicated by arrows in Fig. [j]. In the 
present calculation, however, the van Hove peaks are lo- 
cated at different energies. 

Summarizing this discussion, the tunneling spectra are 
consistent with a sharp van Hove singularity just below 
the Fermi level, and a very broad van Hove singular- 
ity moving to higher binding energy with underdoping. 
Apart from the dip-hump feature, the questions left open 
by the BCS model are the precise energy shift of the 
van Hove singularity of the bonding band, as well as the 
mechanism of its broadening. 



B. Comparison of tunneling and ARPES 

The absence of a sharp van Hove singularity of the B 
band in the tunneling spectra is a clear indication that 
one should go beyond the bare BCS DOS to explain all 
features in the spectra. In principle, it should be possible 
to take spectral functions directly from ARPES data and 
sum them to compare to the tunneling results. However, 
this is complicated by the strong influence of matrix el- 
ements on the photoemission spectra (see e.g. Ref. |3l| ) 
and of the background in the ARPES data. We there- 
fore rely on the phenomenological model of the spectral 
function at low energy described above [with self-energy 
Eq. (jfj)]. This description crudely includes the effects 
of correlation in the marginal Fermi liquid self-energy 
Eo(w). With such an approach we can first directly com- 
pare tunneling and ARPES, and also study the ability of 
the marginal Fermi liquid model to account for tunneling 
data. 

The ARPES data could be fitted using Eq. ([§ with 
«o = 2. In Fig. we compare the calculations for oq = 2 
to experimental tunneling data. Note that the energy 
scale is a factor 2 smaller than in Figs. [I] and |^. As can 
be seen, the calculated curves (dotted lines) do not fit 
the experimental data at all. In particular, the V-shape 
of the spectra at low energy is too narrow. One might 
expect that the excess spectral weight in the gap can be 
suppressed by an anisotropic matrix element. However, 
as shown by the dashed lines, this is not the case. The 
reason is that the low-energy spectral weight is domi- 
nated by the marginal Fermi liquid self-energy, and not 
by the gap function. 

Surprisingly, a good fit can be obtained by an or- 
der of magnitude reduction of «o, to ao = 0.2, as- 
suming an isotropic matrix element. This difference in 
Eq = ImEo fot-ARPES and tunneling suggests, not for 
the first time,t3 that the scattering rate in ARPES is 
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FIG. 2: Tunneling spectra calculated using the phenomeno- 
logical self-energy Eq. (^) and compared to the experimental 
data (circles). Spectra are shown for ao — 2 with the isotropic 
(dotted curve) and anisotropic (dashed curve) matrix element, 
the latter highlighting the (tt, 0) region of the Brillouin zone, 
and for ao = 0.2 with the isotropic matrix element. The 
two curves for ao = 2 are normalized to the spectral weight 
between —300 and +300 meV, while the curve for a = 0.2 
is normalized to the peak height. Ao = (1 + A)~ 1 Ao is the 
renormalized gap value. All spectra are for A = 1 and (3o = 4. 



considerably higher than in tunneling. The explanation 
for this remarkable difference goes beyond our current 
understanding of ARPES and tunneling. 

The van Hove singularity of the B band is shifted closer 
to the Fermi level as a result of band renormalization by 
a factor (1 + A) -1 . Furthermore, the energy dependence 
of Eg leads to a broadening which is more important at 
higher energies. At negative bias this results in a co- 
herence peak followed by the broadened van Hove sin- 
gularity of the B band. However, the energy shift of 
this van Hove singularity as a function of doping is too 
small compa red to the broad feature in the background 
(see Section III A). Of course one can question the va- 
lidity of the-assumption Re Eo = — Xuj at energies above 
~ 50 meVjS meaning that the van Hove singularity of 
the B band is at too low binding energy in our simu- 
lations. Taking this into account would probably lead 
to a better agreement with the broad maximum in the 
background, shifting away from the Fermi level with un- 
derdoping. However, this behavior is inconsistent with 
the sharpness of the dip when the van Hove singular- 
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ity is at much higher binding energy than the dip (as is 
specifically the case for the optimally and underdoped 
samples). 

C. Coupling to a collective mode 

So far, it has not been established whether the cou- 
pling of quasiparticles to the collective mode observed 
in neutron scattering experiments, is sufficient to hawe a 
sizable influence on the spectral functionsEa or not.E2 In 
the following analysis we assume it is, and verify to which 
extent the tunneling spectra are consistent with neutron 
scattering data on the Ok, it) mode. 

In this model [Eq. (|q)], the main effect of the reso- 
nant mode is to enhance the imaginary part £ of the 
self energy b etween the energies Si = — f2 — Ao and 
s-2 = — SI — \AvHs ^r^o> where £ v h s is the band energy 
at the (n, 0) point.EZl A similar effect, although much 
smaller, exists at positive bias. The numerical results 
are displayed in Fig. |^. Focusing first on curves A, which 
represent the DOS of the A band, we notice the sharp co- 
herence peaks. The energy of the latter is a combination 
of the gap value and the van Hove singularity which, as 
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V (mV) 

FIG. 3: Tunneling spectra including interaction with a 
bosonic mode at wave vector (it, it) and energy Q — 5.4 k^T c . 
The coupling constant and correlation length are g = 0.65 eV 
and £ = 2, respectively. The superconducting gap has d-wave 
symmetry and the tunneling matrix element is isotropic. The 
contributions of the A and B bands are shown separately, and 
their sum should be compared to the experimental data. 



in the previous models, lies very close to the gap edges, 
and contributes much to the total spectral weight of the 
coherence peaks. Furthermore, it is crucial for the cre- 
ation of the dip-hump feature, below the coherence peak 
at negative bias. The width of the dip corresponds to the 
energy interval between E\ and £2 where E" is enhanced. 
Below the dip, there is a hump resulting from the scat- 
tering out of the above energy interval. The mode energy 
deferred from neutron scatteringc3cj results in a reason- 
able agreement with the dip in the tunneling spectra (see 
also Ref. |J). However, we caution that its shape is more 
dependent on the precise values of Ao and £ V Hs than on 
the mode energy itself. Furthermore, for the optimally 
and underdoped samples the calculated dip is at slightly 
higher binding energy than in the experimental data, sug- 
gesting that the mode energy is smaller than determined 
from neutron scattering. Finally, the depth of the dip in 
the experimental data indicates an increasing coupling 
constant g with underdoping. 

Let us now turn to the DOS of the B band. Here the 
coherence peaks are small and occur at Ao- The energy 
interval between e\ and £2 is much larger than for the A 
band, and below £2 we see again a hump structure. The 
van Hove peak of the B band now appears inside the dip 
energy range as a very broad maximum. Comparing the 
curves at different dopings, we find that the broad max- 
imum shifts away from the Fermi level with underdop- 
ing. Though this shift is smaller than suggested by ex- 
periment, it is in qualitative agreement with experiment. 
The only disagreement occurs at higher binding energies, 
where the hump coming from the B band is not observed 
in experiment. This hump is probably smeared out by, 
the continuum of spin excitations at higher energies, till 
not taken into account here. 

In Fig. U we compare experimental data to spectra cal- 
culated using different matrix elements. It is clear that 
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FIG. 4: Comparison of the STS spectra including the coupling 
to the (tt, 7r) mode, calculated with different matrix elements. 
Ao = 39 meV, as in Fig. ^| The continued line corresponds to 
the isotropic case Tk = T a , the dotted line to the anisotropic 
case Tfc = T<j, and the dashed line to a partly anisotropic 
matrix element |T fc | 2 = 0.4|T S | 2 + 0.6|T d | 2 . The spectra have 
been normalized to the peak height at —44 meV. 
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the experimental spectrum at low bias shows the V-shape 
typical for a d-wave superconductor, and not the U-shape 
expected for tunneling with a completely anisotropic ma- 
trix element. Though some anisotropy is not to be ex- 
cluded, the tunneling spectra are thus more consistent 
with isotropic than with anisotropic tunneling matrix el- 
ements that suppress all states along the diagonals of the 
Brillouin zone. This conclusion is independent of the ex- 
act model chosen here (it can also be witnessed from, for 
example, the bare DOS in Fig. Q), and only dependent 
on the assumption of d-wave symmetry of the order pa- 
rameter and a small imaginary part of the self-energy at 
low energy. 

IV. SUMMARY 

We have modeled scanning tunneling spectra of 
BSCCO including d-wave BCS superconductivity, band 
dispersions based on recent ARPES data (with bilayer 
splitting), and isotropic as well as anisotropic tunneling 
matrix elements. In addition to this, we have compared 
tunneling spectra to ARPES data via a phenomenolog- 
ical marginal Fermi liquid approach. In a third model, 
we have taken into account coupling of quasiparticles to 
a collective mode with momentum (tt,tv). All numerical 
results have been compared to experimental data over a 
large doping range. 

The simple d-wave BCS model reproduces several gen- 
eral characteristics of the experimental tunneling spec- 
tra, but fails to account for salient features like the dip- 
hump structure and the absence of a sharp van Hove 
singularity from the bonding band. The comparison be- 
tween tunneling and ARPES suggests a much longer life- 
time in tunneling experiments than in ARPES: using the 
marginal Fermi liquid approach, the shape of the coher- 
ence peaks in tunneling spectra is correctly reproduced, 
but only if the lifetime is taken an order of magnitude 
larger than inferred from ARPES data, an observation 
that remains to be explained. This model, specified for 
low energy, does not lead to a satisfactory description of 



the dip-hump structure and the asymmetric background 
as a function of doping. Important improvement is ob- 
tained by including the interaction of quasiparticles with 
a collective mode. With parameters inferred from neu- 
tron scattering experiments, reasonable agreement with 
the tunneling spectra is found, though the energy of the 
mode may be slightly overestimated for the optimally and 
underdoped samples. 

In general, the tunneling spectra are consistent with 
the presence of a sharp van Hove singularity (of the an- 
tibonding band) integrated in the coherence peaks, and 
a broad van Hove singularity (of the bonding band) in 
the background, shifting away from the Fermi level with 
underdoping. These van Hove singularities can be held 
responsible for the asymmetry in tunneling spectra with 
respect to zero bias. The broadening of the bonding band 
van Hove singularity can be largely attributed to the col- 
lective mode. This mode also leads to a dip-hump feature 
in the spectra, which is more pronounced in the under- 
doped than in the overdoped samples. The exact shape 
and depth, however, also depend on van Hove singular- 
ity of the antibonding band and on the superconducting 
gap. 

Finally, the shape of the spectra at low-bias voltage 
is indicative of an isotropic tunneling matrix element. 
We conclude that the tunneling matrix element does not 
have a strong dependence on the (in-plane) wave vector. 
Tunneling spectroscopy therefore probes states along the 
whole Fermi surface, including the diagonals of the Bril- 
louin zone. 
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